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^ ABSTRACT 

>• ' We apply the semi-analytical analysis of the steady nature of line-driven winds 

^ ■ presented in two earlier papers to disk winds driven by the flux distribution of a 

O I standard Shakura & Sunyaev (1973) disk for typical cataclysmic variable (CV) 

Q ■ parameters. We find that the wind critical point tends to be closer to the disk 



in 
o 



surface towards the inner disk regions. Our main conclusion, however, is that 
a line-driven wind, arising from a steady disk flux distribution of a standard 
Ph' Shakura- Sunyaev disk capable of locally supplying the corresponding mass flow, 

O ' is steady. These results confirm the findings of an earlier paper that studied 

' "simple" flux distributions that are more readily analyzable than those presented 

■ here. These results are consistent with the steady velocity nature of outflows 
observationally inferred for both CVs and quasi-stellar objects (QSOs). We find 

^ ' good agreement with the 2.5D CV disk wind models of Pereyra and collaborators. 

■ These results suggest that the likely scenario to account for the wind outflows 
commonly observed in CVs is the line-driven accretion disk wind scenario, as 
suggested early-on by Cordova & Mason (1982). For QSOs, these results show 
that the line-driven accretion disk wind continues to be a promising scenario 
to account for the outflows detected in broad absorption line (BAL) QSOs, as 
suggested early-on by Turnshek (1984), and analyzed in detail by Murray et al. 
(1995). 

Subject headings: accretion, accretion disks — hydrodynamics — novae, cata- 
clysmic variables — QSOs: absorption lines 
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1. Introduction 

As discussed in Pereyra et al. (2004, hereafter Paper I) and Pereyra (2005, hereafter 
Paper II), accretion disks are commonly believed to be present in both cataclysmic variables 
(CVs) and quasi-stellar objects and active galactic nuclei (QSOs/AGNs). A property that 
CVs and QSOs have in common is that both types of objects sometimes present blue-shifted 
absorption troughs in UV resonance lines, giving direct observational evidence for an out- 
flowing wind. Both types of objects also show the existence of a persistent velocity structure 
in their absorption troughs (when present) over significantly long time scales (Paper I). 

In order for a disk wind to account for the wide/broad resonance line absorption struc- 
tures observed in many CVs and QSOs, it must be able to account for the steady velocity 
structure that is observed. The 2.5D time-dependent disk wind models of Pereyra and col- 
laborators, both for CVs and QSOs, have steady disk wind solutions (Pereyra 1997; Pereyra, 
Kallman, & Blondin 1997, 2000; Hillier et al. 2002; Pereyra & Kallman 2003). The ear- 
lier ID disk wind models of Murray and collaborators also find steady disk wind solutions 
(Murray et al. 1995; Murray & Chiang 1996; Chiang & Murray 1996; Murray & Chiang 
1998). 

However, in the literature it has been argued that line-driven disk winds are "intrinsi- 
cally unsteady" because of the increasing gravity along the streamlines at wind base that 
is characteristic of disk winds (see Paper I; Paper II). Since the steady nature of CV and 
QSO wind flows is an observational constraint, whether line-driven disk winds are steady is 
a significant issue. Our objective is to study this issue through a semi- analytical method in- 
dependent of our previous 2.5D numerically intensive simulation efforts (e.g., Pereyra 1997; 
Pereyra, Kallman, & Blondin 2000). 

In Paper I we showed that an increase in gravity at wind base, that is characteristic of 
disk winds, does not imply an unsteady wind solution. We also developed mathematically 
"simple" models that mimic the disk environment and we showed that line-driven disk winds 
can be steady. In Paper II we extended the concepts presented in Paper I, and discussed in 
detail many aspects that were mentioned in Paper I. 

The goal of this paper is to apply the concepts introduced in Paper I and Paper II to 
more realistic models that implement the exact fiux distribution above a standard accretion 
disk. We compare and find good agreement with the steady CV disk wind solutions found 
through numerically intensive calculations by Pereyra (1997) and Pereyra, Kallman, & 
Blondin (2000). 

In §2 we present a general discussion of the reduction of the 2.5D models of Pereyra 
(1997); Pereyra, Kallman, & Blondin (2000) to a series of ID models by neglecting radial 
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terms. We shall see that, although this a rough approximation high above the disk, it is 
quite accurate near the disk surface. We reduce the 2.5D model of Pereyra and collaborators 
to a series of ID equations in §3. In §4 we present the results of our series of ID models for 
different disk radius. We show that the critical point in the inner disk regions is found near 
the disk surface, allowing for an analysis of the solutions with steady wind mass loss rates. 
A summary and conclusions are presented in §5. 

2. General Comments 

The analysis of existence for steady line-driven wind solutions discussed in Paper 1 and 
Paper II correspond to ID models. Thus, to apply them to a 2.5D model, we first reduce 
the 2.5D model to a scries of ID models. To do this we consider the following: the wind in 
CVs tend to be vertical (perpendicular to the disk). There is observational, theoretical, and 
computational evidence for this. In CVs the P-Cygni type line profiles are observed only in 
low inclination systems (close to disk face-on) and with high inferred disk mass accretion rates 
(e.g., Warner 1995). This orientation angle effect was one of the observational properties of 
CVs that led to the early suggestion by Cordova & Mason (1982) that the outfiows detected 
in CVs originate from the accretion disk. 

Prom a theoretical standpoint, if the wind originates from a steady disk, then the radial 
forces at wind base (disk surface) must be at equilibrium, thus the wind must start off 
vertical. We note here that by "vertical" we actually mean "non-radial," since the wind 
starts off with an azimuthal component due to disk rotation. We also note here that although 
the disk wind must in principal start vertical at its base, this does not necessarily imply that 
the wind will continue predominantly vertical. For example, in the case of QSOs/AGNs, 
disk wind models (Murray et al. 1995; Hillier et al. 2002) indicate that the wind tends to 
quickly become radial in nature. 

However, computational 2.5D models of CV disk winds have produced winds that tend 
to be vertical in nature. Theoretical resonance line profiles calculated from these models 
are consistent with observed line profiles in their general form, in the magnitudes of wind 
velocities implied by the absorption components, in the FWHM of the emission components, 
and in the strong dependence with inclination angle (Pereyra 1997; Pereyra, Kallman, & 
Blondin 2000). 

Therefore, near the wind base (disk surface), a CV disk wind at a given radius can be 
well approximated by a vertical ID model. High above the disk rotational forces come into 
play (no longer at equilibrium with the radial component of gravity), streamlines collide. 
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and the density/velocity wind structure will be 2.5D (i.e., 2D with an azimuthal velocity 
component) (Pereyra 1997; Pereyra, Kallman, & Blondin 1997, 2000; Pereyra & Kallman 
2003). Although a ID vertical model is not accurate high above the disk for the above 
reasons, it does give accurate results near wind base. 

As was shown by Castor, Abbott, & Klein (1975) for early type stars, and was discussed 
in Paper II for ID line-driven winds, the wind mass loss rate is determined by the position 
of the critical point. As we show below, for the CV parameters we use in this section [taken 
from the models of Pereyra, Kallman, & Blondin (2000)], at radii where there is significant 
contributions to the wind mass loss, the critical point is near the disk surface. Therefore, 
disk wind characteristics for 2.5D CV models can be addressed. 



The equations of state, mass, and momentum of the 2.5D CV disk wind model of 
Pereyra, Kallman, & Blondin (2000) are, respectively: 



3. 



Reduction to a Series of ID Models 



(1) 



dp , ld{rpVr) , d{pV,) 







(2) 



dt r dr dz 



and 



- 5 - 



dt 



dr 



dVr 
dz 



-P 



GM, 



wd 



dP KeTri'^^z) 



(3) 



+ p X k max 



PHeVth 



dz 



P^ + P^r^ + p-^ + pVz^ = ^ ; 



dt 



dr 



dz 



(4) 



P^T + P^r^- + pVz 



dt 



dr 



dz 



-P 



wd 



(r2 + z2) (^2 + ^2)1/2 



+ p ^ ' ' 



(5) 



Szir.z) , / 
X k max 



PKeVth 



dVz 



dz 



, 10^ 



where P is the pressure, p is the density, b is the isothermal sound speed, Vr, V^, and Vz 
are the corresponding velocity components in cylindrical coordinates, G is the gravitational 
constant, M^^ is the mass of the white dwarf, is the Thomson cross section per mass, c 
is the speed of hght, Vth is the ion thermal velocity, and k and a are the CAK line force 
parameters. J^z and are the corresponding radiation flux components given by: 
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where in turn Q(r) is the rate of energy per area, at a given radius r, radiated by a standard 
disk (Shakura & Sunyaev 1973), that is: 
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where rj is the inner disk radius. In CVs, rj is approximately equal to the white dwarf radius. 
We note that in equations (6)- (9) there is an implicit assumption that rj <^ r/, where r/ is 
the outer disk radius. 

The total disk luminosity C^isk is given by 



^disk = / 4.TTrQ{r)dr , 
J n 

where r/ is the outer disk radius. Thus 
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Considering equations (13) and (14), equation (10) can be rewritten as 
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Additionally, we define two dimensionless functions T(r, z) and (^{r, z): 
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It follows that 
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In reducing the 2.5D model to a series of ID models, for simplicity, we assume an 
isothermal wind. We do not expect this approximation to significantly change the results 
of this work, as is illustrated in the intermediate models of Pereyra, Kallman, & Blondin 
(2000) [see Figures 6 and 7 of Pereyra, Kallman, & Blondin (2000)]. This, on one hand, 
significantly simplifies the analysis, and on the other hand gives a predefined temperature 
structure [namely: b'^{x) — constant] that is a requirement of the ID models of this work 
(see Paper I; Paper II). Thus the equation of state we assume, for the application to the 
2.5D CV disk wind model, is 



(constant) . 



(20) 



We also assume that the streamlines, starting from the disk at a given radius r, are predomi- 
nantly vertical. As discussed in §2, this assumption is justified on observational, theoretical, 
and computational grounds. 

The mass conservation equation is equivalent to the ID models of Pereyra, Kallman, & 
Blondin (1997), namely, for a given r 



-8- 



M^pVAttR^ (constant) , (21) 
where we are determining the existence of ID steady (stationary) solutions, and where 

= r^ + z'^ , (22) 

and M is an integration constant, z is the height above the disk, and r is the radius where a 
given streamhne begins. Thus for each value of r considered, a given ID model is developed. 
We note that in equation (21), M is not the wind mass loss rate (as in Paper I and Paper II), 
but only an integration constant. We have used the symbol M following the notation of 
Paper I and Paper II. However, the wind mass loss rate Mwind can be obtained by (after first 
determining the values of M for different values of r) 

M(r) 

For z <^ r, R"^ ^ {— constant for a given ID model), thus for z <^ r, pV — constant 
[equation (21)] that is the integration of equation (2) assuming stationary solutions and 
neglecting radial velocities. Therefore, equation (21) is accurate near the disk surface and it 
is also accurate towards infinity {z ^ r, ^ z^). However, equation (21) is not accurate 
at intermediate values of z where 2.5D effects become important (Pereyra 1997; Pereyra, 
Kallman, & Blondin 1997, 2000). Using the accurate results near the disk surface, as 
discussed in §2, we study whether wind solutions with steady mass loss rates exist. 

In reducing the 2.5D model to a series of ID models, we additionally neglect the effect 
of a maximum possible value for the line radiation force (Castor, Abbott, & Klein 1975; 
Abbott 1982), that arc represented in the models of Pereyra (1997) and Pereyra, Kallman, 
& Blondin (2000) by including a "max" function in the line force terms [see equation (5)]. 
On one hand, we do not expect this approximation to significantly change the results of 
this work. In the intermediate models of Pereyra, Kallman, & Blondin (2000) this effect 
was found to be most significant in the low density region high above the inner disk region. 
On the other hand, this allows the momentum equation to take a form equivalent to the 
momentum equation of Paper II, that is the form of the momentum equation assumed in 
the ID models of this work. 



Thus, within the assumption of a predominately vertical wind, assuming stationary 
solutions, neglecting radial velocities, and considering equation (5), the momentum equation 
for the ID models (at a given r) becomes 
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where T(r, z) and (^(r, 2;) are given by equations (16) and (17), respectively. The series of ID 
models derived from the 2.5D disk wind model, can be represented within the ID models of 
Paper II, through the following parameterization: 



B(z) = 9M^l-?^rnr,z) . (25) 



A{z) = AttR'^ = 471 {r"^ + z'^) , (26) 
and 7(^) = ^-^^rar,z)k(^Y ' (^7) 

\KeVthJ 



where the body force B{z), the area function A{z), and the line opacity weighted flux 7(2;) 
are defined for a given r. Thus, the mass conservation equation and the momentum equation 
have the form of the equations of Paper II. The isothermal equation of motion is therefore, 

[l- — \AV— = -BA^-tA[^V—] +b^— . 28 
\ V^J dz \M dz J dz ^ ^ 

Defining the characteristic distance Tq, gravitational force Bq, area Aq, and line-opacity 
weighted fiux 70 (see Paper I; Paper II) through 
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the equation of motion becomes 
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where 



x^- ; (31) 
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The parameters s, cu, and m (see Paper I; Paper II) are given by 
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Taking Xq = and q'o = 0, the spatial variable q (see Paper II) 
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becomes 

g = arctan(x) . (38) 
The equation of motion, expressed in terms of q, becomes 
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where the function h{q) is now given by 



h{q) — —ga + Astan{q) 
The nozzle function n for each r is given by (Paper II): 
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Additionally we shall recall the definition of the P function as defined in Paper I and Paper II 



Pico) = 1 - - 
w 



(42) 



that we shall apply in the following section. 



Thus, we have reduced the 2.5D models to a scries of ID models by solving the hydro- 
dynamic equations in the vertical direction and neglecting radial terms. That is, for a given 
radial distance r, wc have a ID model with q = arctan(,T) = arctan(z/r) as the independent 
spatial variable. As discussed above these models are accurate near the disk surface. Since, 
as we see below, the corresponding critical point for inner radii is near the disk surface, the 
ID models for the inner disk regions determine that disk wind flows with steady mass loss 
rates are possible. The models can be used to estimate local mass loss rates (by determining 
local density at the sonic point) as well as velocity and density wind structure near the disk 
surface in the inner region. 



4. Results 

For the ID wind models, we implement the set of parameters used by Pereyra, Kallman, 
& Blondin (2000) corresponding to a CV, namely: 



M^d^O.BMQ ; Cdisk = J^e ; = O.O229i?0 ; 

n^OMRQ ; b = lOkms"^ ; = 2.67 km s"^ 



(43) 



We use the line force parameters also used by Pereyra, Kallman, & Blondin (2000) [and 
Castor, Abbott, & Klein (1975) and Abbott (1982)], namely 



k = 1/3 



a = 0.7 . 



(44) 



As with the analysis of the isothermal CAK stellar wind (Paper 11), we first consider, 
for a given r, the h function for the ID vertical wind model with the flux distribution of a 
standard Shakura & Sunyaev (1973) disk. As discussed in Paper 11, h must be negative 
at the critical point (i.e., h must be negative in order for there to exist values of q, ur, and 
duj/dq that satisfy the critical point conditions). Second, we consider the nozzle function, 
that as we found in Paper II, for the case of an isothermal wind must be a locally increasing 
function at the critical point (i.e., it must have a positive spatial derivative at the critical 
point). Third, for the spatial points where the first two conditions hold, through the nozzle 
function n, the f3 function, the normalized wind mass loss rate m, and the calculation of 
(Puo/dq^ under critical point conditions, as discussed in Paper II, we determine if a local 
solution exists by verifying whether or not the condition 



holds. Fourth, within the set of points that fulfill the first three conditions, we iteratively 
determine the position of the critical point such that the wind reaches sound speed at the 
sonic height of the model (as discussed in Paper II). 

As a consistency check, for a given r, wc superimpose a plot of the nozzle function n 
and the I3{uj)rn function evaluated with the corresponding velocities V{z). As discussed in 
Paper II, in the supersonic wind regime where the h function is negative (Region II), for 
points other than the critical point, the following conditions must hold 



13" {u) fn {Jf + 13' {u) mu" - n"{q) < 



(45) 



/3(a;) 771 < n{q) 



(for u > s and h{q) < and q ^ q^) 



(46) 



and at the critical point 



/3(a;c) m = n{qc) 



(47) 



For r = 2 Tj, Figure 1 shows that the h function is negative from the sonic point to beyond 
1000 times the sonic height. Figure 2 shows that the nozzle function n, corresponding to a 
vertical streamline at r = 2ri, is a monotonically increasing function from the sonic point 
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to beyond 1000 times the sonic height. It is shown, in Figure 3, that equation (45) holds 
from the sonic height to beyond 1000 times the sonic height, showing that local solutions 
to the equation of motion exist throughout the fore mentioned spatial region. By "local 
solution" at a given spatial point, we mean the integration of the equation of motion in the 
vicinity of the given point, assuming the point to be the critical point. Figure 4 shows the 
velocity vs. height obtained upon integrating the equation of motion with the condition of 
achieving the correct sonic height of the model. As a consistency check. Figure 5 verifies 
that equations (46) and (47) hold for the solution obtained. 

Figures 6-10, Figures 11-15, and Figures 16-20 plot the results from an equivalent analy- 
sis to wind streamlines starting at r = 10 r^, r = 20 r^, and r = 50 n respectively. The results 
show that steady line-driven disk wind solutions exist for a standard Shakura-Sunyaev disk. 

A significant result we find is that the critical point tends to be closer to the disk 
surface in the inner disk regions, and farther out from the disk surface in the outer disk 
regions. Physically this is due to the length scales in these different regions. The disk 
wind, at each radii, results from the balancing of the gravitational forces and the radiation 
pressure forces. This balance can be presented quantitatively by the nozzle function defined 
in Paper I and Paper II. The spatial dependence of the nozzle function in turn is crucial in 
determining the critical point position. As larger disk radii are considered, the scale length 
of the corresponding nozzle function also increases, and therefore the critical point position 
increases as well. 

Also, the resulting local wind mass loss rates, for the different radii considered, are 
in good agreement with the 2.5D numerically intensive CV disk wind models of Pereyra 
(1997); Pereyra, Kallman, & Blondin (2000); for example for a radius of lOr^ the local 
wind mass loss rates found in this work is 5.0 x 10~^g s~^ cm^^, and for the same physical 
model parameters the previous numerically intensive 2.5D calculations estimated a value of 
7.1 X 10^'g s^^ cm^^ for the same radius. By confirming the overall earlier results of Pereyra 
and collaborators, we conclude that the likely scenario for the wind outflows in CVs is a 
hne-driven disk wind. This scenario was suggested for CVs early on by Cordova & Mason 
(1982). Pereyra (1997); Pereyra, Kallman, & Blondin (2000) showed that the hne-driven 
disk wind scenario was able to account for the general forms of the C IV AA1549 line proflle 
and its general dependence on viewing angle [for a review on CV properties in general, see 
Warner (1995)]. 

The steady nature of line-driven disk wind that we find here may also be relevant to 

QSO studies. Turnshek (1984) had suggested this scenario early-on based on observational 
constraints. The first serious attempt to study this scenario was by Murray et al. (1995), who 
developed ID streamline models and found that this scenario could account for several QSO 



-14- 



observational features, such as the terminal wind velocities inferred from broad absorption 
lines (BALs) (when present) and the general form of single-trough BALs (when present). As 
discussed in the Introduction, an additional constraint on the BAL region, is that it presents 
a steady velocity structure. Our results here show that the line-driven disk wind scenario is 
consistent with the observed steady velocity structure. 

Thus both observational constraints and theoretical/computational studies to date seem 
to indicate that the line-driven disk wind scenario is a promising one to account for the BALs 
commonly observed in QSOs. 

In future work we plan to study several aspects of this scenario and develop more realistic 
models. 



5. Summary and Conclusions 

In Paper I we had shown that steady wind solutions can exist by using "simple" models 
that mimic the disk environment. These models are more readily analyzable than the more 
detailed models presented here. In Paper II we extended the concepts of Paper I, and dis- 
cussed in detail aspects of the steady/unsteady wind analysis that was presented in Paper I. 
The objective of this work is to determine, in a manner independent of the results of previous 
numerically-intensive 2.5D hydrodynamic simulations, whether steady line-driven disk wind 
solutions exist under the flux distribution of a standard disk (Shakura & Sunyacv 1973). 

Our main conclusion pertaining to the more realistic models presented here is that 
a line-driven wind, arising from a steady disk using the flux distribution of a standard 
Shakura-Sunyacv disk model, is steady. As we had discussed in Paper II, when including 
gas pressures effects, the spatial dependence of the nozzle function continues to play a key 
role in determining the steady /unsteady nature of supersonic line-driven wind solutions. 

These results are consistent with the steady nature of the ID streamline disk wind 
models of Murray and collaborators (Murray ct al. 1995; Murray & Chiang 1996; Chiang 
& Murray 1996; Murray & Chiang 1998). These results also confirm the results of the 2.5D 
steady disk wind models of Pereyra and collaborators (Pereyra 1997; Pereyra, Kallman, & 
Blondin 2000; Hillier et al. 2002; Pereyra & Kallman 2003). 

Another result that we find is that the critical point of the wind tends to be closer to 
the disk surface as smaller radii closer to the inner disk region are considered. 

As we develop more realistic models, we aim towards placing stronger constraints on 
the accretion disk wind scenario, particularly as they apply to QSOs. If the scenario proves 
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tenable, then we may have a well-defined route towards a better understanding of these 
astronomically fundamental objects. 

We wish to thank Kenneth G. Gayley and Norman W. Murray for many useful discus- 
sions. This work is supported by the National Science Foundation under Grant AST-0071193, 
and by the National Aeronautics and Space Administration under Grant ATP03-0104-0144. 
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Fig. 1. — The h function for the ID vertical wind model with the flux distribution of a 
standard Shakura-Sunyaev disk at r = 2ri. A necessary requirement for the critical point is 
that the h function must be negative at that point. This figure shows that all points ranging 
from the sonic point to beyond 1000 times the sonic height fulfill this condition. The physical 
parameters used are: M — 0.6 Mq, Cdisk — ^q: — 0.0229 it!©, = O.Oli?©, h — lOkms"^, 
Vth = 2.67kms~\ A; = 1/3, and a = 0.7 . 
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Fig. 2. — The nozzle function n for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyaev disk at r = 2rj. A necessary requirement for the critical 
point in an isothermal line-driven wind is that the nozzle function must be locally increasing 
(i.e., dn/dq > [dn/dz > 0]) at that point. This figure shows that all points ranging from 
the sonic point to beyond 1000 times the sonic height fulfill this condition. The physical 
parameters used are the same as in Figure 1. 
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Fig. 3. — Local solution existence for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyaev disk at r — 2rf a local solution exists in the vicinity of a 
critical point z if the expression P"m{u;'Y + (3' mu" — n" (solid curve) is negative. This 
figure shows that all points ranging from the sonic point to beyond 1000 times the sonic 
height fulfill this condition. In this plot the terms (3" m{u'Y (dotted curve), P'rhuj" (short 
dash curve), and —n" (long dash curve) are also shown. The physical parameters used are 
the same as in Figure 1. 
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Fig. 4. — Velocity vs. position for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyaev disk at r = 2rj. The critical point position is determined 
with the condition that, upon integration of the equation of motion, the correct sonic point 
position is found. "C" indicates the critical point and "S" indicates the sonic point. The 
physical parameters used are the same as in Figure 1. 
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Fig. 5. — Necessary condition for the global solution existence for the ID vertical wind 
model with flux distribution of a standard Shakura-Sunyaev disk at r = 2 rj for the critical 
point shown in Figure 4: upon the integration of the equation of motion, it must hold that 
Pico) m < n{q) at points other than the critical point [when h{q) < and the wind is 
supersonic], and f3{uj)i7i = n{q) at the critical point. Presented here is the nozzle function 
n (solid curve) and the (3 m function (dotted curve) vs. position. "C" indicates the critical 
point. The physical parameters used are the same as in Figure 1. 
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Fig. 6. — The h function for the ID vertical wind model with the flux distribution of a 
standard Shakura-Sunyaev disk at r = 10 rj. A necessary requirement for the critical point 
is that the h function must be negative at that point. This figure shows that all points 
ranging from the sonic point to beyond 1000 times the sonic height fulfill this condition. 
The physical parameters used are the same as in Figure 1. 
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Fig. 7. — The nozzle function n for the ID vertical wind model with the flux distribution of 
a standard Shakura-Sunyaev disk at r = 10 rj. A necessary requirement for the critical point 
in an isothermal line-driven wind is that the nozzle function must be locally increasing (i.e., 
dn/dq > 0) at that point. This figure shows that all points ranging from the sonic point to 
beyond 1000 times the sonic height fulfill this condition. The physical parameters used are 
the same as in Figure 1. 
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Fig. 8. — Local solution existence for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyaev disk at r = 10 rf. a local solution exists in the vicinity of 
a critical point z if the expression P" m {ui'Y + [3' muj" — n" (solid curve) is negative. This 
figure shows that all points ranging from the sonic point to beyond 1000 times the sonic 
height fulfill this condition. In this plot the terms /3"m(a;')^ (dotted curve), P'rhu;" (short 
dash curve), and —n" (long dash curve) are also shown. The physical parameters used are 
the same as in Figure 1. 
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Fig. 9. — Velocity vs. position for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyaev disk at r = lOrj. The critical point position is determined 
with the condition that, upon integration of the equation of motion, the correct sonic point 
position is found. "C" indicates the critical point and "S" indicates the sonic point. The 
physical parameters used are the same as in Figure 1. 
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Fig. 10. — Necessary condition for the global solution existence for the ID vertical wind 
model with the flux distribution of a standard Shakura-Sunyaev disk at r — 10 ri for the 
critical point shown in Figure 9: upon the integration of the equation of motion, it must 
hold that (3{u!) m < n{q) at points other than the critical point [when h{q) < and the wind 
is supersonic], and /3(a;) m = n{q) at the critical point. Presented here is the nozzle function 
n (solid curve) and the (3 rh function (dotted curve) vs. position. "C" indicates the critical 
point. The physical parameters used are the same as in Figure 1. 
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Fig. 11. — The h function for the ID vertical wind model with the flux distribution of a 
standard Shakura-Sunyaev disk at r = 20 rj. A necessary requirement for the critical point 
is that the h function must be negative at that point. This figure shows that all points 
ranging from the sonic point to beyond 1000 times the sonic height fulfill this condition. 
The physical parameters used are the same as in Figure 1. 
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Fig. 12. — The nozzle function n for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyaev disk at r = 20 r^. A necessary requirement for the critical 
point in an isothermal line-driven wind is that the nozzle function must be locally increasing 
(i.e., dn/dq > 0) at that point. This figure shows that all points ranging from ^ 1.107 times 
the sonic height to beyond 1000 times the sonic height fulfill this condition. However, it is 
also shown that points ranging from the sonic point up to ~ 1.107 times the sonic height do 
not fulfill the condition. The physical parameters used are the same as in Figure 1. 



-28- 




Fig. 13. — Local solution existence for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyaev disk at r — 20ri: a local solution exists in the vicinity of 
a critical point z if the expression P" tti {cu'y + (5' muj" — n" (solid curve) is negative. This 
figure shows that all points ranging from ^ 1.107 times the sonic height to beyond 1000 
times the sonic height fulfill this condition. From ^ 1.107 times the sonic height down to 
the sonic point, the critical point conditions cannot be met (sec Figure 12). In this plot the 
terms l3"rn{u}'Y (dotted curve), P'rhuj" (short dash curve), and —n" (long dash curve) are 
also shown. The physical parameters used are the same as in Figure 1. 
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Fig. 14. — Velocity vs. position for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyaev disk at r = 20rj. The critical point position is determined 
with the condition that, upon integration of the equation of motion, the correct sonic point 
position is found. "C" indicates the critical point and "S" indicates the sonic point. The 
physical parameters used are the same as in Figure 1. 



-30- 





log(z/zJ 



log(z/z,) 




log(z/z,) 




Fig. 15. — Necessary condition for tiie global solution existence for the ID vertical wind 
model with the flux distribution of a standard Shakura-Sunyaev disk at r = 20 for the 
critical point shown in Figure 14: upon the integration of the equation of motion, it must 
hold that Piou) m < n{q) at points other than the critical point [when h{q) < and the wind 
is supersonic], and f3{u!) fn = n{q) at the critical point. Presented here is the nozzle function 
n (solid curve) and the /3 m function (dotted curve) vs. position. "C" indicates the critical 
point. The physical parameters used are the same as in Figure 1. 
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Fig. 16. — The h function for the ID vertical wind model with the flux distribution of a 
standard Shakura-Sunyaev disk at r = 50 rj. A necessary requirement for the critical point 
is that the h function must be negative at that point. This figure shows that all points 
ranging from the sonic point to beyond 1000 times the sonic height fulfill this condition. 
The physical parameters used are the same as in Figure 1. 
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Fig. 17. — The nozzle function n for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyaev disk at r = 50 rj. A necessary requirement for the critical 
point in an isothermal line-driven wind is that the nozzle function must be locally increasing 
(i.e., dn/dq > 0) at that point. This figure shows that all points ranging from 2.11 times 
the sonic height to beyond 1000 times the sonic height fulfill this condition. However, it is 
also shown that points ranging from the sonic point up to ~ 2.11 times the sonic height do 
not fulfill the condition. The physical parameters used are the same as in Figure 1. 



-33- 



log(z/z,) 




log(z/z,) 




1 2 3 4 1 2 3 

log(z/z,) log(z/zJ 



Fig. 18. — Local solution existence for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyacv disk at r = 50 rj: a local solution exists in the vicinity of 
a critical point z if the expression ji" m (c<j')^ + 13' muo" — n" (solid curve) is negative. This 
figure shows that all points ranging from ^2.11 times the sonic height to beyond 1000 times 
the sonic height fulfill this condition. From fa 2.11 times the sonic height down to the sonic 
point, the critical point conditions cannot be met (see Figure 17). In this plot the terms 
P"m{u;'y (dotted curve), P'rhu;" (short dash curve), and —n" (long dash curve) are also 
shown. The physical parameters used are the same as in Figure 1. 
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Fig. 19. — Velocity vs. position for the ID vertical wind model with the flux distribution 
of a standard Shakura-Sunyaev disk at r = 50rj. The critical point position is determined 
with the condition that, upon integration of the equation of motion, the correct sonic point 
position is found. "C" indicates the critical point and "S" indicates the sonic point. The 
physical parameters used are the same as in Figure 1. 
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Fig. 20. — Necessary condition for the global solution existence for the ID vertical wind 
model with the flux distribution of a standard Shakura-Sunyaev disk at r = 50 for the 
critical point shown in Figure 19: upon the integration of the equation of motion, it must 
hold that l3{uj) m < n{q) at points other than the critical point [when h{q) < and the wind 
is supersonic], and l3{uj) m = n{q) at the critical point. Presented here is the nozzle function 
n (solid curve) and the /? m function (dotted curve) vs. position. "C" indicates the critical 
point. The physical parameters used are the same as in Figure 1. 



